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We present a systematic method to introduce free parameters in sets of mutually unbiased bases. 

In particular, we demonstrate that any set of m real mutually unbiased bases in dimension N > 2 
admits the introduction of (m — l)IV/2 free parameters which cannot be absorbed by a global 
unitary operation. As consequence, there are m = fc+ 1 mutually unbiased bases in every dimension 
N = k 2 with k 3 /2 free parameters, where k is even. We construct the maximal set of triplets 
of mutually unbiased bases for two-qubits systems and triplets, quadruplets and quintuplets of 
mutually unbiased bases with free parameters for three-qubits systems. Furthermore, we study the 
richness of the entanglement structure of such bases and we provide the quantum circuits required to 
implement all these bases with free parameters in the laboratory. Finally, we find the upper bound 
for the maximal number of real and complex mutually unbiased bases existing in every dimension. 

This proof is simple, short and it considers basic matrix algebra. 
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I. INTRODUCTION 

Mutually unbiased bases (MUB) have an ubiquitous 
role in quantum mechanics. They are useful to generate 
quantum key distribution protocols m, detection of en¬ 
tanglement 3|, quantum random access codes 0, dense 
coding, teleportation, entanglement swapping and covari¬ 
ant cloning (see [B] and references therein). Furthermore, 
a maximal set of MUB allows us to univocally reconstruct 
quantum states [7i|. On the other hand, entropic certainty 
00 and uncertainty relations m, have been considered 
for MUB. Such important applications have motivated an 
enormous effort to understand the underlying structure 
behind incomplete mm and complete mm sets of 
MUB. In particular, incomplete sets of MUB have an 
important role in Bell inequalities |14j . uncertainty rela¬ 
tions [15] and locking of classical correlations in quantum 
states fTSHTgj . Despite of the important advance done 
for complete sets of MUB in prime [7] and prime power 
[T3] dimensions, incomplete sets of MUB seem to be much 
more challenging. Indeed, the full classification of all pos¬ 
sible sets of MUB for 2-qubit systems has been recently 
done m and the multipartite case is poorly understood. 
The lack of a deeper understanding of mutually unbiased 
bases seems to be the absence of a suitable mathemat¬ 
ical tool. Indeed, for three or more qubits systems it 
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is not known the existence of quadruplet of MUB hav¬ 
ing free parameters and a few triplets were accidentally 
found [50]. In this work, we enlighten this area of research 
by presenting a systematic method to introduce free pa¬ 
rameters in incomplete sets of MUB. As consequence, we 
demonstrate that any set of m real MUB existing in any 
dimension N admits the introduction of free parameters 
with our method. We also show that our construction 
is not restricted to the consideration of real bases. In¬ 
deed, we illustrate our method by explicitly constructing 
the maximal set of triplets for 2-qubits systems and sev¬ 
eral triplets, quadruplets and quintuplets of MUB having 
free parameters for 3-qubits systems. All of these cases 
involve complex MUB. Furthermore, we analyze the en¬ 
tanglement structure of such sets of MUB and provide 
the quantum circuit required to implement all these sets 
in the laboratory. 

This work is organized as follows: In Section II we 
present a short introduction to mutually unbiased bases, 
the link to complex Hadamard matrices and we resume 
the state of the art of MUB with free parameters. In Sec¬ 
tion III we present our method to introduce free param¬ 
eters in incomplete sets of MUB. In Section IV we prove 
that any set of real MUB admits the introduction of the 
maximal number of parameters allowed by our method. 
In Section V we construct triplets, quadruplets and quin¬ 
tuplets of MUB having free parameters for three qubit 
systems and study the entanglement structure of each 
case. In Section VI we resume our main results, conclude 
and discuss some open questions. Additionally, we illus¬ 
trate our method by explicitly solving the simplest case 
of triplets of MUB for 2-qubits systems (see Apendix|A|. 
The explicit construction of a quadruplet and a quintu¬ 
plet of MUB having free parameters for three qubit sys- 
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terns is provided in Appendix[B] In Appendix[C]we derive 
the quantum circuit required to generate every quadru¬ 
plet and quintuplet of MUB presented in this work. Fi¬ 
nally, in Appendix [D] we find a simple and short proof 
for the upper bound of the maximal number of real and 
complex MUB in every dimension by considering basic 
matrix algebra. 


II. MUTUALLY UNBIASED BASES AND 
COMPLEX HADAMARD MATRICES 

In this section, we present some fundamental proper¬ 
ties about mutually unbiased bases (MUB) and complex 
Hadamard matrices (CHM) required to understand the 
rest of the work. For a complete review about MUB 
and CHM we suggest [21] and [22] . respectively. Two or¬ 
thonormal bases {0y}y=o,...,JV—l and {4’k}k=o,...,N-i de¬ 
fined in C N are mutually unbiased if 

\{<t>j,^k)\ 2 = (1) 

for every j, k = 0,..., N — 1. In general, a set of m > 
2 orthonormal bases are MUB if every pair of bases of 
the set is MUB. A set of m MUB is called extensible if 
there exist an m + 1th basis which is mutually unbiased 
with respect to the rest of the bases. It has been shown 
that m = N + 1 MUB exist for N prime [7] and prime 
power Ha. Additionally, maximal sets of MUB can be 
constructed in prime power dimensions by considering 
gaussian sums and finite fields For any other 

dimension N = p r ip r 2 ■ ■ -p r k k ( p\ 1 < p r £ < • • ■ < p r k ) the 
maximal value for m is not known and the lower bound 
to > pY is provided by the maximal number of fully 
separable (i.e., tensor product) MUB [2B]. Additionally, 
in dimensions of the form N = k 2 it is possible to find 
to = k + 1 real MUB by considering orthogonal Latin 
squares m- Let us arrange the bases and {tf>k} 
in columns of unitary matrices B i and B 2 , respectively. 
Thus, if Bi and B 2 are MUB we have 

B\B 2 = H, (2) 

where H is a CHM. An N x N matrix H is called a 
complex Hadamard matrix (CHM) if it is unitary and 
all its complex entries have the same amplitude 1/y/N. 
For example, the Fourier matrix (F)v)jfc = -A=e 2,r *- 7 ' fc / JV 

is a CHM for every N, where i = y/—l. Two CHM Hi 
and H 2 are equivalent if there exists permutation matri¬ 
ces Pi, P 2 and diagonal unitary matrices Di, D 2 such that 
H 2 = PiDiHiD 2 P 2 . Therefore, MUB and CHM are close 
related: any set of to MUB Si = {Hi,..., B m } is unitary 
equivalent to a set S 2 = {I, Hi, ..., H m - 1 }, where I rep¬ 
resents the computational basis and Hi,... ,H m _ 1 are 
CHM. Indeed, the unitary transformation that connects 


Si with S 2 is B\. That is, 

B\(Si) = {B\Bi,b\b 2 , ..., B\B m } 

= {I ,Hi,...,H m _i} 

= S 2 , (3) 

where we used Eq.|2j). Alternatively, H) ; (Si) also pro¬ 
vides an analogous result for k = 2,..., to. The full 
classification of CHM and MUB has been solved up to 
dimension N = 5 (see 23 and [19], respectively). For 
N = 6 both problems remain open despite a consider¬ 
able effort made during the last 20 years H2E2H3ZI- The 
problems also remain open for any dimension N > 6. For 
example, they are open in the prime dimension N = 7, 
where a single 1-parametric family of complex Hadamard 
matrices is known [38] and a maximal set of 8 MUB is 
known [7] but incomplete sets of MUB are not yet char¬ 
acterized. Indeed, it is still open the question whether a 
triplet of MUB having free parameters exist in dimension 
N = 7. 

Let us summarize the state of the art about sets 
of MUB having free parameters. First, any set of 
to MUB in prime dimension N = p of the form 
{I, F p , Ci ,..., C m - 2 } is isolated, where F p is the Fourier 
matrix and {Ci,..., C m - 2 } are circulant CHM [39]. A 
complex Hadamard matrix is isolated if there is no fam¬ 
ily of CHM connected with it 22- By family we un¬ 
derstand a set of inequivalent complex Hadamard ma¬ 
trices depending on some free parameters. We extend 
the same definition to sets of MUB: a set of to MUB is 
isolated if there is no family of to MUB connected with 
it. For example, any set of to < A^ + 1 MUB in di¬ 
mension N = 2,3 and N = 5 is isolated. In dimension 
N = 4 there is a unique 3-parametric triplet of MUB of 
the form {I, F^\x), H(y, z)} and quadrupets and quin¬ 
tuplets of MUB are isolated [T9] . In dimension N = 6 a, 
1-parametric triplet of MUB exists [99] . Moreover, two- 
parametric triplets of the form {I, Fg (x, y),H(x, y)} ex¬ 
ist for any x,y € [0,2ir) and seem to be unextendible for 
any pair x, y Also, in dimensions N = 9 [401 and 

N = 4fc 1-parametric triplets of MUB can be defined by 
considering cyclic n-roots pfO] . 

All the above sets of MUB with free parameters were 
found by taking advantage of special properties holding 
in specific dimensions. The existence of quadruplets of 
MUB having free parameters is still unknown in every 
dimension, as far as we know. In the next section we 
present a systematic method to introduce free parameters 
in sets of to MUB in dimension N. 


III. MUB WITH FREE PARAMETERS 

A set of r > N vectors {v k } C C N has associ¬ 
ated a Gram matrix G £ C rxr , where Gij = ( Vi,Vj), 
i,j = 0, ...,r — 1 and Rank(G) = N. Reciprocally, 
from the Gram matrix G of size r and rank N we can 
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always find a set of vectors {v' k } such that = {v'^v'j) 
and v' k £ C N . The set of vectors {ufc} and {u{}, associ¬ 
ated to the same G, are connected by means of a unitary 
transformation. Thus, they define the same geometri¬ 
cal structure in the complex projective space CP^' 1 . 
The vectors {?;{} can be found from G by considering 
the Cholesky decomposition , i.e., to find the unique up¬ 
per triangular matrix L having positive diagonal entries 
such that G = L^L. Thus, the r vectors {v' k } are given 
by the r columns of L, where we only have to consider 
entries of the first N rows of L (the rest are null because 
of the rank restriction). In this work, we are particularly 
interested to study Gram matrices associated to a set of 

That is, 
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( 4 ) 

Note that this matrix naturally defines a structure of m 2 
square blocks of size N, each of them defined by a uni¬ 
tary matrix of the form HjHj , where i,j = 0,..., m — 1 
and Hq = I. The Cholesky decomposition of this Gram 
matrix G is given by 


L = 


(l Hi H 2 
On On 0w 


\ Ojv 0 N 0 n 


H m -1 \ 
On 


On J 


( 5 ) 


where On are zero matrices of size N. So, the set of m 
MUB is clearly given by {I, Hi, H 2 , ..., H m -i} which cor¬ 
responds to the first block of rows of G (see Eq. 0 )- This 
important property substantially simplifies our method. 
In a previous work, we found the most general way to in¬ 
troduce free parameters in pairs of columns (or rows) of 
any complex Hadamard matrix in every dimension [i'l] . 
A free parameter can be introduced in two columns C\ 
and C 2 of a CHM if and only if C\ o C 2 £ R N , which only 
holds for N even. Here, the circle denotes the (entrywise) 
Hadamard product, that is, (Gi o C 2 )j = ( Ci)j(C 2 )j, 
j = 0, ...,N — 1. Pairs of columns (or rows) satisfy¬ 
ing this property were called equivalent to real pairs (ER 
pairs). The construction of the CHM having free param¬ 
eters is very simple: 

CONSTRUCTION III.l (gU). Given an ER pair of 
columns {Ci,C 2 } we introduce a free phase e la in the 
jth entries ( Ci)j and ( C 2 )j if (C* o C 2 )j < 0 for j = 
" .V 1. 


Here, the asterisk denotes complex conjugation. Note 
that ° C 2 )j is the inner product between the 

column vectors G 4 and C 2 , which has to be zero by def¬ 
inition of CHM. Let us exemplify this method by intro¬ 


ducing two free parameter in the Fourier matrix 
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( 7 ) 


Here, we considered the ER pairs of columns {Ci,C 3 } 
and {C 2 ,C 4 } to introduce the paramaters a and ft, re¬ 
spectively. Note that a and p are aligned in the same 
rows. A set of N/2 ER pairs producing aligned free pa¬ 
rameters in matrices of size N are called aligned ER pairs. 
The remarkable property of aligned ER pairs is that one 
of the parameters is always linearly dependent (if we con¬ 
sider the equivalence of CHM defined above). Thus, the 
parameter p is linearly dependent and the Fourier family 
has only one relevant parameter [22] . 

Let us now extend this method to the construction of 
MUB with free parameters. Here, the key ingredient is 
the generalization of the concept of ER pairs: a set of 
two columns {Gx,G 2 } of a Gram matrix of m MUB in 
dimension N is called a generalized ER pair (GER) if C 4 o 
C 2 £ U. mN . The following result, natural generalization 
of Construction |III.1| is the main result of this work: 

PROPOSITION III.l. Let G be the Gram matrix of 
a set of m MUB in dimension N and suppose that it has 
A f GER pair of columns, where both vectors of each GER 
pair belong to the same block of columns. Then, the set 
of MUB admits the introduction of Af free parameters. 


Proof. For simplicity let us first consider the case of 3 
MUB (to = 3) in dimension N, where the Gram matrix 
is given by 


G= Hi 


H 4 

I 


Hi H^Hx 


H 2 

h\h 2 

II 


( 8 ) 


and suppose that G has a GER pair of columns {Gi, Cj} 
in the same block of columns (i.e., Int[*/AT] = Intg/IV], 
where Int means integer part). Therefore, N of the prod¬ 
ucts ( Ci)k(Cj)k are zero because of the corresponding 
identity block I and only 2 N values of these products play 
a role in Gi o Cj. Thus, a free parameter can be intro¬ 
duced in both columns Gi and C 3 by applying Construc¬ 
tion fTTTTTl to the 2 N dimensional subvectors of C,- and Go- 


having 27V non-zero entries. Note that after introducing 
the parameter the hermiticity of G is destroyed. In order 
to restore it we have to apply the same method to the 
GER pair of rows {/?,;, Rj} which always exists because 
of the hermiticity of G. This lead us to a 1-paranxetric 
set of matrices satisfying i) G(a ) = U(a)GU\a) and 
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ii) \G(a)ij\ = \Gij\ for any a £ [0, 2n). Note that i) 
holds for any U whereas ii) is strongly dependent on our 
construction. Here, U(a) and W(a) represent the in¬ 
troduction of a free parameter in columns two and rows, 
respectively. Furthermore, G(a) and G(0) have the same 
eigenvalues for any a £ [0,27r) so G(a ) is a 1-parametric 
set of Gram matrices defining a 1-parametric set of to 
MUB in dimension TV. If G has J\f GER pairs then we 
can introduce A f free parameters in the same way. The 
generalization to any to > 3 is straightforward from the 
above explanation. □ 


Let us emphasize the importance of considering both 
vectors of a GER pair in the same block of columns: sup¬ 
pose that we choose a GER pair formed by columns of 
different blocks (e.g. the 5th and 9th columns of the 
Gram matrix given in the example Eq.(A4|) and we in¬ 
troduce a free parameter. Despite of this action generates 
a genuine Gram matrix the set of 3 bases would be not 
composed by MUB with free parameters. This is sim¬ 
ple to understand because of the free parameter would 
appear in a single vector of the second and third basis. 


In the particular case of to = 2 our Proposition III.l 


is reduced to Construction |III.T] which has been derived 
in a previous work [41] , That is, to introduce free pa¬ 
rameters in a pair of MUB is equivalent to introduce free 
parameters in a CHM, as suggested by Eq.([3]). It worths 
to mention that GER pairs of columns belonging to the 
first block of G (i.e., {Ci,Cj} with i,j < TV) always pro¬ 
duce parameters that can be absorbed in global unitary 
transformations. Roughly speaking, in this case the pa¬ 
rameters do not appear in the inner product of vectors 
of two different bases (see Eq.|4)). The explicit construc¬ 
tion of the maximal set of triplets of MUB for 2-qubits 
systems is given in Appendix[Aj Also, a triplet, a quadru¬ 
plet and a quintuplet of MUB having free parameters for 
3-qubits systems are given in Appendix [B] We encourage 
to the reader to have a close look to the examples in order 
to clearly understand our method. 


IV. FAMILIES STEMMING FROM REAL MUB 

As we have shown, Proposition [TTTTT] allows us to intro¬ 
duce free parameters in Gram matrices of MUB having 
GER pairs. In this section, we demonstrate that every 
set of to of real MUB in dimension TV > 2 allows the in¬ 
troduction of the maximal number of parameters allowed 
by GER pairs: 

PROPOSITION IV.1. Any set of m real MUB in di¬ 
mension TV > 2 admit the introduction of Nm/2 free pa¬ 
rameters. Furthermore, (to— l)TV/2 of these parameters 
cannot he absorbed by global unitary transformations and, 
at most, one of them is linearly dependent. 

Proof. Every pair of columns belonging to the same block 
is clearly a GER pair. Therefore, there are Nm/2 GER 
pairs allowing the introduction of Nm/2 free parameters. 


The rest of the proof is stra ghtforward (already explained 
in the proof of Prop. III.l). □ 


Furthermore, note that there are many different ways 
to define the GER pairs and so many different families of 
MUB can be constructed. Precisely, there are (^) differ¬ 
ent ways to define GER pairs in each of the to— 1 blocks of 
columns (the first block only provides unitary equivalent 
MUB). That is, a total of (to — 1)TV(TV — l)/2 different 
ways. We do not know how many ways are inequiva¬ 
lent for TV > 4. As we mentioned before, in dimensions 
N = k 2 it is possible to construct m = k + 1 real MUB. 
By combining this result with the above proposition we 
have the following result: 


COROLLARY IV. 1 . In every dimension N = k 2 there 
exists to = k + 1 MUB admitting k 3 /2 free parameters. 


Here, we consider k > 1 and thus kN parameters can¬ 
not be absorbed by a global unitary transformation. In 
dimension TV = 4 there exists to = 3 real MUB and thus 
we can introduce k 3 /2 = 4 free independent parameters, 
where the GER pairs are aligned and thus one of the 
four parameters is linearly dependent (see Appendix |A|. 
Here, the 12 possible ways to introduce free parameters 
produce equivalent sets m- Such triplets are also equiv¬ 
alent to the solution found in Appendix [A] Sets of real 
MUB are not the only cases where a maximal number of 
parameters can be introduced with our method. Indeed, 
in the next section we construct sets of MUB with free 
parameters in the 3-qubit space by considering complex 
MUB. 


V. MUB FOR THE 3-QUBIT SPACE 

In dimension TV = 2 3 there is a maximal number of 9 
MUB. Indeed, four maximal sets having a different en¬ 
tanglement structure have been constructed [23]: 

Si = (2,3,4) S 2 = (1,6,2) 

S 3 = (0,9,0) <S 4 = (3,0,6) (9) 

where the first, second and third entries denote the num¬ 
ber of fully separable, biseparable and maximally entan¬ 
gled bases, respectively. These sets of MUB are not equiv¬ 
alent under Clifford operations but they are equivalent 
under general unitary transformations. In this section 
we focus on the construction of pairs, triplets, quadru¬ 
plets and quintuplets of MUB having free parameters and 
stemming from elements of <S 4 . We have chosen this par¬ 
ticular set of MUB because it contains the highest num¬ 
ber of maximally entangled bases and, consequently, it 
has potentially important applications in quantum infor¬ 
mation theory. 


From considering Prop. |III.l| we find the following results 
for every subset of m. < 9 MUB of <S 4 : 
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#MUB(m) 

H Param. 

Example 

2 

3 

(I, Hi}* 

2 

4 

(I, HJ 

3 

7 

{I ,HcH 2 y 

3 

8 

{I,H 1 ,H 2 } 

4 

4 

{i,h 1 ,h 2 ,h 3 } 

5 

0 


5 

4 

{I,H 1 ,H 2 ,H 3 ,H 5 y 

6-9 

0 

{I, H,, H :l , H k , Hi, H m } 


TABLE I: MUB with free parameters in dimension N = 8. 
The asterisk means that GER pairs are aligned, which pro¬ 
duces one linearly dependent parameter. Note that a set of 
MUB in general allows many different choices for GER pairs, 
and some of them produce non-alligned GER pairs. In Ap¬ 
pendix [B] we detailedly explain all possible sets of m MUB 
having free parameters. 


(i) Every pair of MUB {I, Hf} C S 4 admits the intro¬ 
duction of 4 free parameters, where i = 1,..., 8 . 

(ii) Every triplet of MUB {I, Hj, H 3 ] C 1 S 4 admits the 
introduction of 8 free parameters for every i 7 ^ j = 
1 ) ■ • • 8 - 

(in) Some quadruplets of MUB {I, Hj, Hk, Hi} C £4 ad¬ 
mit the introduction of 4 free parameters in only 
one of the bases, whereas the rest of the quadru¬ 
plets do not admit free parameters. 

(iv) Every quadruplet admitting 4 free parameters can 
be extended to a quintuplet of MUB having 4 free 
parameters. The extension of quadruplets to quin¬ 
tuplets is not unique. 

(v) Every set of 6 < m < 9 MUB do not admit free 
parameters. 

The maximal number of parameters that can be intro¬ 
duced for every m is provided in Table [B]of Appendix [B] 
where we present the proof of the above results. 

The free parameters of all quadruplets and quintuplets 
of MUB described in (Hi) and (iv) can be generated in 
the laboratory by considering 7 different quantum circuits 
(56 cases, 8 cases per circuit; see Appendix [C|), which in¬ 
volve local and Toffoli gates. The generation of the fixed 
set of bases, i.e. the set £> 4 , requires different a circuit 
[42] which involves local, non-local controlled-phase and 
Toffoli gates [45]. Therefore, the sets of MUB with free 
parameters are generated by a composition of two dif¬ 
ferent quantum circuits. The explicit expression of the 
quantum circuits is provided in Appendix [C] 

The entanglement structure of the sets of MUB pre¬ 
sented in Appendix[C]is very interesting. For example, let 
us consider the quintuplet {I, Hi (a), H 2 , H 3 , H 5 } C <S 4 
(explicitly constructed in Appendix H- For simplicity, 
let us assume that the 4 parameters aq to aq are iden¬ 
tical (a). Here, Hi(0) is a maximally entangled basis, 



PA 


PB 

PC 


tt/4 

a 


*/2 


FIG. 1: Purity of the reductions to Alice (pa), Bob (pb) and 
Charlie (pc) for all the states of the basis Hi(a). As we can 
see, Bob and Charlie are as entangled as they can for every 
a whereas Alice is maximally entangled with Bob-Charlie for 
a = 0 (GHZ state) and separated for a = k/2, where Bob 
and Charlie share a Bell state. 


in the sense that every vector of the basis is equiva¬ 
lent (i.e., up to local unitary operations) to the GHZ 
state | GHZ) = (|000) + 1111))/x/2. On the other hand, 
Hi ( 7 r/ 2 ) is a biseparable basis. Indeed, every vector of the 
basis is equivalent to \cj>) = |0)(|00) + |ll))/v / 2. That is, 
Alice is separated and Bob and Charlie share a maximally 
entangled Bell state. For any 0 <a< 7 r/ 2 we have an 
intermediate amount of entanglement between Alice and 
Bob-Charlie whereas Bob and Charlie are as entangled as 
possible for any a. That is, two of the three parties (Bob 
and Charlie) saturate the maximal amount of entangle¬ 
ment allowed by the monogamy of entanglement of three- 
partite systems. Indeed, for any value of the parameter 
a the single qubit reductions pb and pc are maximally 
mixed. It is highly non-trivial the fact that the basis H i(0) 
(maximally entangled) and Hi(tt/ 2) (fully separable) can 
be continuosly connected without loosing the unbiasity of 
the quintuplet of MUB for any value a £ [ 0 , 7 r/ 2 ). In Ap¬ 
pendix [C] we show that, in this case, the parameter a is 
fully controlled by local unitary operations generated by 
Bob (see quantum circuit G). This means that Bob has 
full control on the entanglement existing between Alice 
and Bob-Charlie when we are restricted to keep the unbi¬ 
asity of the quintuplet. Table [B] in Appendix [C] shows all 
possible entanglement structures that can be find from 
quadruplets and quintuplets with free parameters arising 
from £> 4 . The purity of the reductions pa, Pb and pc as a 
function of the free parameter a for the above quintuplet 
is depicted in Fig. 0 


VI. SUMMARY, CONCLUSION AND 
DISCUSSION 

We have presented a systematic way to introduce free 
parameters in sets of m mutually unbiased bases in di- 
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III. 1). In particular, for 


mension N (see Proposition 
TO = 2 our method is reduced to introduce free param¬ 
eters in complex Hadamard matrices (see Construction 


III. 1 and also our previous work EH)- We proved that 


any set of m real mutually unbiased bases existing in 
any dimension N > 2 admit the introduction of free 
parameters. Furthermore, in every dimension N = k 2 
there are k + 1 mutually unbiased with kN/2 free pa¬ 
rameters, where k is even. We have found the maximal 
set of triplets of mutually unbiased bases with free pa¬ 
rameters for two qubit-systems (see Appendix |A|) . Also, 
we constructed pairs, triplets, quadruplets and quintu¬ 
ples of mutually unbiased bases having free parameters 
for 3-qubit systems (see Appendix [0. Such sets are con¬ 
structed from subsets of the maximal set of complex MUB 
S 4 (see Eq.0). That is, our construction is not restricted 
to sets of real MUB. Additionally, we provided the com¬ 
plete set of quantum circuits required to implement all 
such quadruplets and quintuplets (see Appendix JC|. Fi¬ 
nally, we presented a new proof for the upper bound of 
the maximal number of real and complex mutually unbi¬ 
ased bases existing in every dimension. This short and 
simple proof only involves basic algebra (see Appendix 

0 - 


The analysis provided in Section [V] for 3-qubits can be 
easily extended to a higher number of qubits. In order to 
do this we have to consider the construction of maximal 
sets of mutually unbiased arising from Galois fields [25] , 
Such construction is a generalization of the set £ 4 , where 
every basis is a real Hadamard matrix multiplied by a di¬ 
agonal unitary matrix containing 4th roots of the unity. 
In such cases we can always define GER pairs (see Prop. 


of m MUB. Finally, let us present some open issues: ( i ) 
Find the subset of triplets, quadruplets and quintuplets 
of MUB considered in Appendix [C] such that they are 
extendible to 9 MUB, and ( ii ) Is it possible to construct 
maximal sets of MUB with free parameters in some di¬ 
mension? This question seems to have a negative answer 
for every N. However, a formal proof is only known for 
every N < 5 m\- 


III. 1) and, therefore, introduce free parameters in subsets 
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Appendix A: MUB in dimension four 


Let us consider the simplest case where our method 
can be applied. The construction of 771 = 2 MUB hav¬ 
ing free parameters is reduced to find a family of CHM 
and, thus, here we consider 777 = 3. In dimensions N = 2 
and N = 3 complex Hadamard matrices are isolated and, 
consequently, any set of MUB in such dimensions is iso¬ 
lated. On the other hand, in dimension 4 there is a fam¬ 
ily of CHM (see Eq. 0). So, the first case corresponds to 
N = 4. A fixed triplet of MUB for N = 4 is given by: 


ill 


( l 0 0 o\ 
0 10 0 
0 0 10 
V 0 0 0 1 / 




(11 1 1 

11-1-1 
1-1 7-7 

\ 1 — 1 — i i ) 


H , = 


/ 1 1 1 1 \ 
1 11 - 1-1 
y /2 -11 1-1 

V 1 -1 1 - 1 / 


(Al) 


(A2) 


(A3) 


The Gram matrix G associated to this set is given by 1 /2 
of the following matrix 


/ 2 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 \ 

0 

2 

0 

0 

1 

1 

-1 

-1 

1 

1 

-1 

-1 

0 

0 

2 

0 

1 

-1 

7 

— 7 

-1 

1 

1 

-1 

0 

0 

0 

2 

1 

-1 

— 7 

7 

1 

-1 

1 

-1 

1 

1 

1 

1 

2 

0 

0 

0 

1 

1 

1 

-1 

1 

1 

-1 

-1 

0 

2 

0 

0 

1 

1 

-1 

1 

1 

-1 

— 7 

7 

0 

0 

2 

0 

7 

—7 

1 

1 

1 

-1 

7 

—i 

0 

0 

0 

2 

—7 

7 

1 

1 

1 

1 

-1 

1 

1 

1 

— 7 

7 

2 

0 

0 

0 

1 

1 

1 

-1 

1 

1 

7 

—7 

0 

2 

0 

0 

1 

-1 

1 

1 

1 

-1 

1 

1 

0 

0 

2 

0 

\1 

-1 

-1 

-1 

-1 

1 

1 

1 

0 

0 

0 

2 ) 


(A4) 


Here, we have 6 GER pairs of columns and rows given by 
{l-2;3-4;5-6;7-8;9-10; 11-12}. Note that the perfect match 
between GER pairs of columns and rows is given to the 
fact that G is hermitian. The introduction of free param¬ 
eters into the GER pairs of columns {Ci, C^}; {C3, C4} 
and rows {i?i, R 2 }; {-R 3 , -R 4 } implies that and 

HIH 2 do not depend on the parameters (See Eq.0). 
As consequence, these MUB are unitary equivalent for 
any value of the two parameters. Thereby, they are not 
interesting for us. In general, we do not introduce free 
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parameters in the first N columns (and rows) of Gram 
matrices of m MUB in dimension N as we already ex¬ 
plained in Section |lll[ From considering the remaining 4 
GER pairs above defined we easily generate the following 
MUB with free parameters: 


H 2 (a,/3) 


H 3 (l, 5) 



< 1 

1 

1 

1 

1 

1 

1 

-1 

-1 

7^ 

e ia - 

e ia 

ie l P 

—ie 1 ^ 


\e ia - 

e ia 

-ie 

ie l P 


( 1 

1 

1 

1 \ 

1 

1 

1 

-1 

-1 

71 

- e G 

e 

e iS 

_ e iS 


V e i7 

—e* 7 

e iS 

-e u ) 


(A5) 


(A6) 


Note that our method can be considerably simplified by 
introducing the free parameters in a reduced region of 
the Gram matrix (A4). This is because G contains much 
more information than the set of 3 MUB (see Eq. [5]). 
Precisely, we can restrict our attention to introduce free 
parameters in the first 4 rows according to the existing 
pairs of GER and Proposition |III.1| The rest of the rows 
give us the explicit expression of the inner products be¬ 
tween the elements of the different bases. In general, we 
can restrict our attention to introduce parameters in the 
first N rows of G when we consider m MUB in dimension 
N. Note that this implies to only consider the GER pairs 
of the last (m — 1)IV columns (without considering GER 
pairs of rows). 

The f our parameters a,/3, 7 , 5 appearing in Eq.(A5l 
and Eq.(A 6 l cannot be absorbed by global unitary op¬ 
erations. However, they are aligned, so one of them can 
be absorbed in a global phase of a vector of the canonical 
basis Hi. Therefore, we find the following 3-parametric 
continuos triplet of MUB in dimension N = 4: 


{I,H 2 (a,0),H 3 ( 7 ,6)}. (A7) 

This has been reported as the most general triplet of 
MUB that can be constructed in dimension N = 4. 
Quadruplets and quintuplets of MUB do not allow free 
parameters in dimension 4 and, consequently, they are 
isolated [T5] . 


Appendix B: MUB in dimension 8 

In this appendix we construct the maximal number of 
triplets, quadruplets and quintuplets of MUB having free 


parameters from £4 (see Eq.Q). The key result is pro¬ 
vided in Table [B] where we present the complete set of 
GER pairs for 54 . 


How to read Table [B] 

(i) The cell associated to column Hk and row T/j 
contains all the GER pairs allowed by the triplet 
{I ,Hj,H k }. 


(ii) The notation i-j-k-l means that every possible com¬ 
bination of 2 non-repeated indices determine a 
GER pairs; that is, {Ci,Cj}, { Ci,C k }, {Ci,Cj}, 
{Cj,Ck}, {Cj,Ct} and (G\., C/} are GER pairs. 


(in) The semicolon (;) separates complementary sets of 
GER pairs (i.e., for {i-j-k-l;fi-u-K-ij} mixtures of 
Graco-Latin indices do not form GER pairs). 


To construct quadruplets or quintuplets of MUB we have 
to find the intersection of the sets of GER allowed by all 
the subsets of triplets. If there is no intersection then 
free parameters cannot be introduced. Let us construct 
a triplet of MUB: 

Suppose we want to introduce free parameters in the 
triplet {H\, H 2 ,H 3 } (see Eq.| 8 ])). In order to introduce 
free parameters in H 1 we have to find common GER pairs 
in the cells associated to H\Hi (i.e., column 2, row 3 of 
Table [bJ {l-2-5-8;3-4-6-7}) and h\h 1 (i.e., column 2, 
row 4: {l-2-4-7;3-5-6-8}). This is equivalent to find GER 
pairs appearing in the Gram matrix of {Hi, H 2 , H 3 }. 
Thus, the unique set of common GER pairs are given by 
{Ci,C 2 }, {C 3 ,C 6 }, {C 4 ,C 7 } and {C 5 ,C 8 }. Analogously, 
we can find the GER pairs for the second and third block 
of G; that is, {C U C B }, {C 2 ,C 5 }, {C 3 ,C 4 }, {C 6 ,C 7 } 
and {Ci, C 4 }, {C 2 ,C 6 }, {C 3 ,C 8 }, {C 5 ,C 7 }, respectively. 
Thus, we are in conditions to introduce 12 free param¬ 
eters in the Gram matrix of the fixed set {Hi, H 2 , H 3 }. 
As noted in Appendix [Aj the Cholesky decomposition al¬ 
lows us to simplify the introduction of free parameters 
by only considering the first N rows of G. Thus, our 12 
parametric set of m = 3 MUB is given by 













h 2 

Hg 

Hi 

Hg 

Hg 

H 7 

H s 

Hi 

- 

l-3-4-8;2-5-6-7 

l-4-5-7;2-3-6-8 

l-2-4-7;3-5-6-8 

l-4-5-6;2-3-7-8 

l-6-7-8;2-3-4-5 

l-2-3-7;4-5-6-8 

l-2-5-8;3-4-6-7 

Hi 

l-2-5-8;3-4-6-7 

- 

l-3-4-8;2-5-6-7 

l-2-3-6;4-5-7-8 

l-2-5-8;3-4-6-7 

l-2-5-8;3-4-6-7 

l-4-5-7;2-3-6-8 

l-2-4-6;3-5-7-8 

H\ 

l-2-4-7;3-5-6-8 

l-6-7-8;2-3-4-5 

- 

l-4-5-6;2-3-7-8 

l-6-7-8;2-3-4-5 

l-3-4-8;2-5-6-7 

l-2-4-8;3-5-6-7 

l-3-4-8;2-5-6-7 

H\ 

l-3-5-7;2-4-6-8 

l-4-5-7;2-3-6-8 

l-2-4-6;3-5-7-8 

- 

l-2-3-6;4-5-7-8 

l-4-5-6;2-3-7-8 

l-3-5-8;2-4-6-7 

l-4-5-7;2-3-6-8 

Hi 

l-2-3-6;4-5-7-8 

l-2-4-6;3-5-7-8 

l-2-5-8;3-4-6-7 

l-3-5-7;2-4-6-8 

- 

l-3-5-7;2-4-6-8 

l-6-7-8;2-3-4-5 

l-6-7-8;2-3-4-5 

Hi 

l-3-4-8;2-5-6-7 

l-2-3-7;4-5-6-8 

l-3-5-6;2-4-7-8 

l-2-5-8;3-4-6-7 

l-3-4-8;2-5-6-7 

- 

l-2-5-6;3-4-7-8 

l-3-5-6;2-4-7-8 

H\ 

l-6-7-8;2-3-4-5 

l-2-5-8;3-4-6-7 

l-6-7-8;2-3-4-5 

l-6-7-8;2-3-4-5 

l-2-4-7;3-5-6-8 

l-2-3-6;4-5-7-8 

- 

l-2-3-7;4-5-6-8 

Hi 

l-4-5-6;2-3-7-8 

l-3-5-6;2-4-7-8 

l-2-3-7;4-5-6-8 

l-3-4-8;2-5-6-7 

l-3-5-7;2-4-6-8 

l-2-4-7;3-5-6-8 

l-3-4-6;2-5-7-8 

- 


TABLE II: GER pairs required to construct any subset of m MUB steeming from the maximal set of 9 MUB S 4 . 


and 


f 


a 2 , «3,or) — 


V 


H 2 (Pl, /? 2 , 03 1 Pi) — 


^ 3 ( 71 , 72 , 73 , 74 ) = 


—1 

—l 

1 

1 

—* 

1 

l 

— l 

—i 

—i 

—i 

i 

i 

—i 
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l 

-ie iai 

ie iai 

ie iOL2 - 

-ie ia3 - 

ie iot4 

-ie ia2 

ie ia * 

ie ia4 

ie iai 

—ie lai 

ie ia2 

ie ia3 - 

ie 1014 

—ie 1012 

— ie ta3 ie iai 

e iai 

_ e i°n 

-e ia 2 

—e iot3 - 

-e iQ4 

e *“ 2 

e ia 3 

gia 4 

e iai 

-e iai 

e ia2 

-e ia3 

1 

_g^ 4 e *Q :2 glOL 3 

—e iot4 

-1 

-1 

1 

-1 

1 

1 

-1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

/ -1 

1 

-1 

-1 

1 

1 

1 

-i \ 


—ie 1 ^ 2 

—ie 1 ^ 3 

ie 1 ^ 3 

ie 

ie 1 ^ 4 

—ie 1 ^ 4 

ie l @ 1 

—i 

i 

i 

i 

i 

—i 

—i 

—i 

—e i/Sl 

_ e i02 

gifts 

_ e i/3s 

e if>2 

_gifti 

e i/3 4 

gi /? 1 

—e i/?1 

e i02 

e i03 

_ e i/3s 

— e 4 ^ 2 

e i/3 4 

_gi/3 4 

gift 

—i 

—i 

i 

i 

—i 

i 

i 

—i 

ie 1 ^ 1 

—ie 1 ^ 2 

ie I/3s 

—ie 1 ^ 3 

ie 1 ^ 2 

ie 

—ie 1 ^ 4 - 

-ie I ^ 1 

V 1 
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1 

1 

1 

1 

1 

1 / 


/ ie^ 1 

ie 112 

ie 113 

—ie* 71 

- ie 174 

—ie 112 

ie 174 

—ie 113 \ 

-e i71 

e il2 

e* 73 

e* 71 

e* 74 

— e* 72 

—e il4 

— e* 73 

-1 

1 

-1 

-1 

1 

1 

1 

-1 

—i 

—i 

i 

—i 

i 

—i 

i 

i 

—e* 71 

—e il2 

gi73 

e ni 

—e il4 

e il2 

e* 74 

— e^ 73 

-ie i71 

ie il2 

—ie il3 

ie i71 

—ie il4 

- ie i72 

ie il4 

ie il3 

—i 

i 

i 

— * 

— i 

i 

—i 

i 

V 1 

1 

1 

1 

1 

1 

1 

1 / 


(Bl) 


(B2) 


(B3) 


Here, only 8 parameters are relevant because 4 of them 
do not appear in the inner products, and thus they can 
be absorbed by a global rotation (i.e., or a’s or /3’s or 7 ’s 
can be considered as zero without loosing of generality). 
Moreover, given that the parameters are aligned we have 
7 independent parameters. This triplet can be straight¬ 
forwardly extended to a quadruplet of MUB by adding 


the computational basis Hg. In order to construct a quin¬ 
tuplet we have to find a suitable extra basis. One way 
to do this is by considering Hg. For this choice we have 
a non-empty set of GER and 4 parameters can be intro¬ 
duced in Hi (see Table [B| column starting with Hi). The 
remaining four bases I,H 2 ,H 3 and Hg are fixed, where 
H 2 = H 2 (0, 0,0,0), H 3 = H 3 ( 0,0, 0,0) and Hg is given 
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by 


1 

-1 

1 

1 

-1 

-1 

-1 

1 

1 

-1 

-1 

1 

1 

1 

-1 

-1 

—i 

—i 

—i 

i 

i 

—i 

i 

i 

—i 

—i 

i 

i 

—i 

i 

i 

—i 

i 

—i 

—i 

—i 

—i 

i 

i 

i 

—i 

i 

—i 

i 

—i 

i 

—i 

i 

-1 

-1 

1 

-1 

1 

1 

-1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


(B4) 



We encourage to the reader to verify that there is no in¬ 
tersection between the set o GER pairs for Hg , H 3 and 
H 5 (see Table E and, consequently, no more free pa¬ 
rameters can be introduced. Therefore, the 4 parametric 
quintuplet is given by 

{I, H\(a\, 0% aq), H 2 ,Hz, H$}. (B5) 


Circ. A 

Circ. B 

Circ. C 

Circ. D 

Circ. E 

Circ. F 

Circ. G 

1246 

1234 

1236 

1345 

1245 

1237 

1235 

1257 

1278 

1258 

1367 

1267 

1248 

1268 

1347 

1368 

1348 

1468 

1478 

1358 

1378 

1356 

1467 

1456 

1578 

1568 

1457 

1567 

2478 

2358 

2378 

2346 

2347 

2368 

2348 

2568 

2457 

2567 

2357 

2356 

2467 

2456 

3468 

3456 

3467 

2458 

3458 

3567 

3457 

3578 

5678 

4578 

2678 

3678 

3568 

4678 


TABLE III: Quantum circuit required to construct every 
quadruplet and quintuplet of MUB steeming from £ 4 . The 
4 numbers ijkl denote the quintuplet {I, Hi, Hj, Ht, Hi} 
whereas the red-bold number (online version) denotes which 
basis carries the parameters. Quadruplets are constructed by 
removing any basis from quintuplets. 


FIG. 2: Quantum circuit required to construct the fixed set 
of 9 MUB S 4 [42]. 

CHM Hi, Hj, Hk and Hi. The fifth basis is Hg = I which 
is implicit in the table. Thus, the quintuplet associated to 
ijkl is given by {I, Hi, Hj, H^, Hi}. From removing any 
basis we get a quadruplet having free parameters. The 7 
quantum circuits are given by 

QUANTUM CIRCUIT A: 


R(a) 


QUANTUM CIRCUIT B: 


- R(a) 

Appendix C: Quantum circuits for quadruplets and 

quintuplets of MUB - 


Every quadruplet and quintuplet of MUB with free pa¬ 
rameters in dimension 8 shown in this work was generated 
from considering subsets of S 4 (see Eq.Q). Such sets can 
be implemented in the laboratory by considering suitable 
quantum circuits. First, we have to consider the gener¬ 
ation of the fixed bases S 4 and then the introduction of 
the parameters. Therefore, the full quantum circuit is 
the composition of two circuits. The generation of the 
set <S 4 is given by the quantum circuit depicted in Fig¬ 
ure [ 2 ] This circuit was recently derived [32]. The free 
parameters can be introduced by considering 7 quantum 
circuits (A to G). Table [B] shows the quantum circuit 
required to generate every quadruplet and quintuplet of 
MUB. Here, every set of four numbers ijkl denote the 4 


QUANTUM CIRCUIT C: 


R(a) 


QUANTUM CIRCUIT D: 
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QUANTUM CIRCUIT E: 


♦ 





parameters appear according the following table: 


Circuit 


B C D E 


G 


(C2) 

Rows 1234 1256 1357 1467 2367 2457 3456 
Here, every set of 4 numbers determine the rows where 
the free parameters are introduced in the bases. For ex¬ 
ample, the quintuplet {I, H 3 ,H 2l H 3 , H 5 } provided in Ap¬ 
pendix [B] allows the introduction of parameters in H\. 
The fixed set of 5 bases is generated by the quantum cir¬ 
cuit given in Fig(2] According Table |B| the free parame¬ 
ters can be introduced in Hi by considering the quantum 
circuit G. The aligned parameters in H i appear in the 
rows 3,4,5 and 6 (according Table C2|. We encourage 
readers to verify these properties from the explicit ex¬ 
pressions of the quintuplet {I, Hi, H 2 , H 3 , H 5 } provided 
in Appendix [B| 


Appendix D: Maximal number of MUB: a simple 
proof for the upper bound 


Here, we present an independent proof for the upper 
bound of the maximal number of MUB in real and com¬ 
plex Hilbert spaces. 

PROPOSITION D.l. In dimension N there are m,R < 
N/ 2 + 1 and me < N + 1 MUB for real and complex 
Hilbert spaces, respectively. 

Proof. Let G be the Gram matrix of to complex MUB in 
dimension N. Then, 


QUANTUM CIRCUIT G: 


Go G^ = ^1 + 1, 


(Dl) 


♦ 



R(a) 


where 


R(a) = 




(Cl) 


and a £ [0,27r). Every quantum circuit (from A to G) 
introduces aligned free parameters in one of the basis of 
every quintuplet. The selection of the quadruplet or quin¬ 
tuplet has to be according Table [B] as we already ex¬ 
plained and exemplified in Appendix [5] The aligned free 


where the mN x mN matrix 1 has to diagonal blocks of 
size N consisting by the null matrix and the non-diagonal 
blocks of size N are equal to the unit matrix 1 (i.e., every 
entry of 1 is 1). The matrix I of Eq.(Dl I is the identity 
matrix of size mN. From matrix theory it is known that 
given A, B > 0 we have Rank(AoR) < Rank(A)Rank(I3). 
From consi derin g A = B^ = G(n,d), the above inequal¬ 
ity and Eq.(Dl I we have mN — (to — 1) < N 2 or, equiva¬ 
lently, to < N + 1. For the real case, we have the tighter 
inequality Rank(A o A) < |Rank(A)[Rank(A) + 1] [45] . 
From combining this inequality with Eq.(Dl) we have 
m.R < N/2 + 1. □ 


This proof was inspired in the derivation of the upper 
bound of the maximal number of vectors in Equiangular 
Tight Frames [44] . 
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